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Abstract 

We consider the time evolution of TV bosonic particles interacting 
via a mean field Coulomb potential. Suppose the initial state is a 
product wavefunction. We show that at any finite time the correlation 
functions factorize in the limit N ^ oo. Furthermore, the limiting one 
particle density matrix satisfies the nonlinear Hartree equation. The 
key ingredients are the uniqueness of the BBGKY hierarchy for the 
correlation functions and a new apriori estimate for the many-body 
Schrodinger equations. 
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1 Introduction 

It is well-known that the Schrodinger equation governs the dynamics of non- 
relativistic quantum systems. The basic existence and uniqueness of the 
Schrodinger equation were studied extensively for many systems, including 
the important Coulomb systems. Despite these efforts, there are few cases 
that the solutions to the many-body Schrodinger equation can be described 
in reasonable details. One such case is a system of N weakly interacting 
bosons, or more precisely, bosons interacting with a pair potential whose 
coupling constant is of order N~^. 

For such a system, we expect that the potential acting on any single 
particle is approximately generated by the average density of all other par- 
ticles. Therefore, there are no correlations generated by the dynamics and 
all correlations of the state at a later time are almost all due to the initial 
correlations. In the simplest scenario, i.e., when the initial state has no cor- 
relation, the correlations at a later time should be negligible as well. We 
consider initial states of N bosons that are described by products of one- 
particle wave functions. Therefore, our goal is to show that, if the initial 
A^-body wave function is a product of a one-particle wave function, then it is 
well-approximated by a product of some one-particle wave function at later 
times as well. Assuming this is correct, it is easy to see that the one particle 
wave function satisfies a nonlinear Schrodinger equation whose nonlinearity 
is given by the convolution of the potential with the mass density of the one 
particle wave function, i.e., V * |V'p. This equation is typically called the 
Hartree equation. In the special case when the pair interaction is Coulomb, 
it is called the Schrodinger-Poisson system. As usual for the Schrodinger 
equation, we can always cast it in the form of an equation for the density 
matrix. We shall use the same name for both setups. 

This problem has a long history and has been considered by many au- 
thors. For regular pair-potentials, it was solved by Hepp [5] in the context 
of field operators and by Spohn [9] using density matrices. Hcpp's method 
requires differentiability of the potential, while Spohn uses only the bound- 
edness of the potential. Ginibre and Velo [3] have greatly extended these 
results to include singular potentials (including Coulomb). However they 
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worked in the second quantized framework and the quasi-free initial state 
used in [3] describes special classes of excitations around the product state. 
In particular, this method does not apply to wave-functions with a fixed 
number of particles. 

The basic tool to analyze the many-body dynamics is the BBGKY hierar- 
chy. We notice that the time derivative of the fc-particle correlation function 
(more precisely, /c-particle density matrix) of an A'^-body state is given by 
the k + l particle correlation function for each k = 1,2, . . . N — 1. This gives 
a hierarchy of coupled evolution equation for all k correlation functions up 
to N, called the Schrddinger hierarchy. Without further limiting procedure, 
this system is just a rewriting of the original Schrodingcr equation. If we 
take the formal N ^ oo limit, it converges to an infinite system of equa- 
tion, called the BBGKY hierarchy, for the fc-particle density matrices for 
k = 1,2, . . . . For a product initial state this infinite hierarchy has a trivial 
solution of product form built up from the solution of the Hartree equation. 
Although the A^-body hierarchy has a unique solution, the infinite hierarchy 
may, in principle, admit more than one solution. Thus our task is to establish 
the convergence of the Schrodinger hierarchy to the BBGKY hierarchy as 
N ^ oo, and prove the uniqueness of the solution to the BBGKY hierarchy. 

Both the BBGKY and the Schrodinger hierarchy can be put into systems 
of integral equations. For bounded potential, these systems can be solved 
by iteration which converges in the trace norm of the density matrices. Us- 
ing this idea on the Schrodinger hierarchy, Spohn proved that the fc-point 
correlation functions of the N body system are given by convergent power 
series with an error estimate uniformly in N. Thus one can take the limit 
N oo and one obtains that the fc-point correlation functions are indeed 
given by the products of the solution to the nonlinear Hartree equation. 

Bardos et. al. have followed a different route [1]. They showed that any 
t/;*-limit point of the solutions to the A^-body hierarchy satisfies the infinite 
BBGKY hierarchy. This proof requires the potential to be bounded below 
and in -|-L°°. If, in addition, the potential is bounded, then the sequence 
of correlation functions is shown to be 7i;*-precompact in the trace norm. 
The uniqueness of the solution to the infinite hierarchy is established by 
controlling the trace norm similarly to [9]. In particular, this result estab- 
lishes the convergence to the BBGKY hierarchy for the repulsive Coulomb 
systems which was not covered in Spohn's work. 

The main question is the uniqueness of the BBGKY hierarchy in the 

case of singular potentials. Furthermore, the convergence in the attractive 
Coulomb case should be resolved since it describes the important gravita- 
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tional system. Technically, all the uniqueness methods rely heavily on the 
boundedness of the potential via the estimate Tr|y7| < ||F||ooTr|7| for a den- 
sity matrix 7. The Coulomb potential, i.e., V{x) = ±|x|~^, is unbounded, 
but Hardy's inequality allows us to control it in i^^-norm, i.e., 

Tr \Vj\ < C{ Tr V7V + Tr 7 ) . 

Therefore the solution to the BBGKY hierarchy should be unique under 
the right Sobolev norm. Since we have correlation functions of arbitrary 
number of particles, the Sobolev space we choose is the iterative Sobolev 
space (see Section 3 for the precise definition) which is weaker than the usual 
one but sufficient for the uniqueness. This poses the problem of estimating 
the iterative Sobolev norm for correlation functions. This estimate can only 
be obtained from the original A^-body Schrodinger equation for which we 
know only the conservation of mass and energy. Notice that in [9] and [1] 
the only estimates obtained are trace norm and norm bounds, which 
are consequences of the conservation of mass and energy for the A/"-body 
Schrodinger equation. 

Since we can only establish the uniqueness of the BBGKY hierarchy in 
the iterative Sobolev norm, we need to establish an apriori estimate that the 
A;-correlation functions are bounded in such a norm. It should be emphasized 
that to establish such estimates directly for the correlation functions of a 
Coulomb system is very difficult due to the \x\~^ singularity. However, it is 
feasible for a cutoff Coulomb systems with an dependent cutoff. The main 
idea here is to take the energy to a higher power and deduce from there the 
estimate on the iterative Sobolev norm. Notice that the third power of the 
potential is already not in L^. Hence it does not even define a meaningful 
operator in the usual sense. This is part of the technical reason we need to 
perform a truncation for the Coulomb singularity. 

We now need to control the original evolution and the cutoff dynamics. 
Here we use the conservation of the norm to control the difference of these 
two evolutions. Since to control an iV-body system in general produces a 
factor N, our cutoff has to be rather small. The various restrictions on the 
cutoff scale finally give the choice of order o(iV~^/^). 

This work was partly inspired by the work of [1] and partly by the ques- 
tion posed by J. Yngvason regarding the derivation of the time dependent 
Gross-Pitaevskii equation from the many-body Schrodinger equation. Un- 
fortunately, our method, as it stands, still cannot prove the convergence 
to the Gross-Pitaevskii equation. Another possible motivation for studying 
high density bosons with Coulomb interaction is that electrically charged 
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ions may have bosonic statistics. 

This work was supported by NSF grants DMS-9970323 (L.E.) and DMS- 
0072098 (H.-T. Y.) The discussions with C. Bardos, F. Golse, N. Mauser 
and J. Yngvason are gratefully acknowledged. The results of this work and 
[1] were jointly outlined in [2]. 

2 Many-body Schrodinger evolution with mean field 
potential 

Consider a system of A^-bosons weakly interacting with a Coulomb potential. 
The dynamics of such a system is governed by the Schrodinger equation with 
the Hamiltonian 

1 ^ 1 

e=i i<j 

defined on ^f^i L^(R'^), where V{x) = ±fi\x\~^ with > 0. The wave 
functions are symmetric functions of A'' variables. Since the Hamiltonian is 
symmetric, the wave function at the time t will be symmetric as long as the 
initial data is symmetric. 

The Schrodinger equation is given by 

idt^N,t = HN'^N,t, (2.1) 

with the initial data specified at the time t = 0. The equation can be 
solved expHcitly by "^N,t = e~'*^^^Ar,o- Let jN,t = ''"^jvt be the projection 
operator in L^(R^) associated with the wave function '^N,t- The Schrodinger 
equation is equivalent to the operator equation 

idtlN,t = [HN,lN,t] ■ (2.2) 

This is a more general setup since it allows a general density matrix not 
coming from a wave function. The n-point density matrix of 7jv,t is defined 
by 

..X„;x;,...<):=/W.„...x„,..„....x.) (2.3) 

X*jv,f(,Ti, ...x'„, Xn+l, ■ ■ ■ XN)dXn+l ■ ■ ■ dXN 

if n < N and 7]^^ := otherwise. The normahzation is (n < N) 

TV 7^ = 1 . (2.4) 
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It is a simple calculation that the n-point density matrices of the solution 
to the Schrodinger equation satisfy the finite hierarchy 



(2.5) 



k „i 

i^N,k{*)'lNoi^i^---^k;x[,...x'k) + {-i)Y] / dsUN,k{t-s) 

e=i 

J dxk+i(v{xi - Xk+i) - V{x'i - Xfe+i)^ 



N -k 
N 



x1n^^\xi, ...Xk, Xk+i;x[, ...x'k, Xk+i) 



(k) ■ (fc) 

where UN,k{t)l = e~^*^N -ye^^^N and 



k k 



If we take the limit N ^ oo and neglect all lower order terms, this system 
converges to the BBGKY hierarchy given by the following infinite system of 
equations for density matrices 

j^''\xi,...Xk;x[,...x'k) (2.6) 

k ..I 

= i^k{tho\xi,...Xk;x[,...x'k) + {-i)y2 dsUk{t-s) 
J (v{xe - xk+i) - V{x'^ - xk+i)^ 
X 7?^"^^ {xi,...Xk,Xk+i;x'i,...x'i,;,Xk+i )dxk+i 



/here W,(t)7 := [UU^~''^'^^ 



7 [Uj=ie-''^'/^), = 1,2,.... This 
is the integral form of an infinite system of evolution equations with initial 
data To := (7^^^ 7^^\ . . . ) . 



Suppose the initial data is of product form and given by 

N 

*7V,o(a;i, . . . xn) ■■= JJ i^oixj) 



(2.7) 
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for some initial onc-particle wave function ^pQ. We always assume that ipo 
is normalized in L^(R^), i.e. H^olb = 1- For this initial data, the /c-point 
density matrix, is, simply the tensor product 7,^5 of the one-point 
density matrix 7^^, which is the projection matrix onto V'o- One can check 
that if the initial data to the BBGKY hierarchy is of product form, then 
there is a special solution to the BBGKY hierarchy which is of product 
form and the one-particle density matrix satisfies the nonlinear Schrodinger 
equation 



If we denote the kernel by 74(0;, x'), then (2.8) is equivalent to 



idtjt{x,x') 



1 



2 

+ / dz 



A,. - A,./ 



lt{x,x') 



(2.8) 



(2.9) 



V{\x - z\) - V{x - z) -ft{z,z)'yt{x,x') . 



We can put this equation into a more familiar form. If 7o(a;, x') := '4'oix)'4'oix'), 
then 7((x,x') := tpt{x)ijjt{x') where satisfies 

idtMx) = —^xM^) + {Jv{x- z)\Mz)\^dz^Mx) (2.10) 

with initial data '0t=o = V'o- This equation was studied extensively by Gini- 
bre and Velo. In particular, the equation preserves the norm and the 
energy and we have [4] 



sup \\ipt\\H^ < 00 
t>0 



(2.11) 



if the norm of the initial condition ||V'o||_ffi is finite. Although the equa- 
tion is nonlinear, the normalization condition HV'olb = 1 can be assumed 
without loss of generality at the expense of changing 

Therefore, if we can justify the limiting procedure and prove the unique- 
ness of the BBGKY hierarchy, the evolution of the weakly interacting A'^- 
bosons can be understood by a one-body nonlinear Schrodinger equation. 
We first describe the topology for the limiting procedure. 

Denote by H the Hilbert space L^(R^) and let jC{H) be the set of 
bounded operators with operator norm || • ||. We let C^{H) be the set of 
trace class operators with the norm ||7||i := Tr|7|. The set of density ma- 
trices, C^, is defined as the subset of nonnegative self-adjoint trace class 
operators. Let /C := 1C{H) be the set of compact operators equipped with 
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the operator norm. It is well-known (Theorem VI. 26 in Vol.1, of [7]) that the 
dual space of the compact operators is the space of trace class operators, i.e., 
(/C, II • ID* = (C^, II • 111). This gives rise to the variational characterization of 
the trace norm: 



\A\\i = sup 

K&K.{H) : \\K\\=1 



TV AK 



(2.12) 



The u'*-topology on is induced by a family of seminorms A — > |Tr AK\ 
that are indexed by the family of compact operators K e K. 

Denote by H'^'^ := (g)*Li = 0^=1 -^^(R^) the fc-tensor product of 
^•^(R'^). We can define the trace norm on H^^ and extend the duality from 
A; = 1 to all k. Define the space 

CX) 

C:={r = (7^,7(2),...) : j('^')ejC\H^'')] = l[£\H^''). 

k=l 

We equip this set with the product of the it;*-topologies on each component. 
Thus the convergence in C is characterized by the following property: r„ — > T 
as n — ^ 00, if for each k wc have 7!^^ — > 7*-'^'-' in w* sense in Jl^{H'^^). 

Let (t) := (1 + i^)V2. We define the set 

00 

L-(R+,(i)-Mt,C) ■.= \{L^{K^,{t)-^dt,C\H®^)) 

k=l 

i.e., the set of functions 7 W(t) : R+ ^ C^{H'^^^^) with supt(t)-i||7(^)(*) 111 < 
00. We equip this set with the product of the ty*-topologies on each factor. 
On each factor the tf;*-topology is given by seminorms 



^W(i)^ / Tr K{t)-i^^\t) 
Jo L 



'0 

where K{t) G L1(R+, {t)dt,IC{H'^^)), i.e. \\K{t)\\{t)dt < 00. 
Let 

For any one-particle wave function i/j we define 

n 

7i/, (^1) ■ ■ ■ ^nj Xi, . . . Xj^) := 7^(xj, Xj) . 
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Suppose ipt G //^(R'^) is a solution to the nonlinear Schrodinger equation 
(2.10), then the collection of density matrices 

is a solution to (2.6). We now state the main theorem: 

Theorem 2.1. Let V{x) = ±iJ,\x\~^ be the repulsive or attractive Coulomb 
potential with some /x > 0. Assume that V'o G H^(R^) and let ipt be the 
solution of (2.10). Let Fjv,* be the solution to (2.5) with initial condition 
^Nfl '■= r^o- -^s N ^ oo, we have 

in the w* topology o/L°^(R_|_, {t)^^dt,C). In other words, for each k 

(k) (k) 

in the weak* topology o/L°°(R+, {t)-^dt, jC^iH^'')). 



Remark 1: Note that our theorem requires the initial one-particle wave func- 
tion to be in although the natural space for the nonlinear Schrodinger 
equation (2.10) is H^. 

Remark 2: Our proof works for more general potentials as well, for example 
one can consider V{x) = yl(a;)|a;|~^ -|- B{x) with A,B ^ S(R?) Schwarz 
class. The control on high derivatives is necessary because of the commutator 
estimate (5.13). 



We close this section by an important remark. The equations of the 
BBGKY hierarchy are written for the kernels of the density matrices, how- 
ever we interpret them as equations for density matrices even though we 
sometimes use the traditional kernel notation. This point of view allows us 
to circumvent the problem that apriori the kernels axe only functions de- 
fined almost everywhere, hence setting x'j^^^ = x^j^i in the interaction term 
requires an extra argument. In this paper we interpret the integral 

j V{xe - Xk+i)j^'''^^^ {xi,...Xk, Xk+i;x[, ...x'^;, Xk+i)dxk+i 

as the partial trace 

where V is the multiplication operator by V{xi—Xk+i). We will always verify 
that V^^'^'^^^ is a trace class operator on L^(R^(*^"'"^)), hence the partial trace 
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is well defined. In the Appendix we give the precise definition of the partial 
trace and we collect a few basic information. 

Convention: Throughout the paper the letter C will refer to various con- 
stants depending only on /x. 

3 Sobolev spaces of density matrices 

We recall that C{H), C^{H) and IC{H) denote the set of bounded, trace 
class and compact operators on the Hilbert space H, respectively. Also, 
£} C C}{H) denote the set of density matrices. We have the standard 
inequality 



and the variational characterization of the trace norm (2.12). We also recall 
that \A + B\ < |^| + |5| is noitrue in general, butTr|^ + S| < Tr|A| + Tr|S| 
is true, following from (2.12), (3.1). Moreover, if we have two Hilbert spaces, 
Hi, H2, then the following inequality is valid for partial traces 



In the Appendix we give the proof of this inequality and we also collect a 
few basic facts about partial traces. 

We define the analogue of the Sobolev norm for trace class operators 
as follows: 

||r||^i := Tr \STS\ 

where S := {I — A)^/^. We introduce the space of operators = H^^H) 
as: 



and let := H} n C^. Notice that ||7||^i = Tr S-fS = Tr7 + Tr(-A)7 for 



We define a weak* topology on and for this purpose we identify this 
space with a dual space. We let 



AB\\i < \\A\\ \\B\\i 



(3.1) 



IVi IV2 ^ < IVi,2 \A\ . 



(3.2) 





and we equip this space with the norm 
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Lemma 3.1. With the notations above 

(Aii-iur = (HMi-ii«i). 

The proof of this Lemma is found in the Appendix and it is similar to the 

standard proof of (/C, || -11)* = {C^, \\ • ||i) outhncd in [7]. Using this duahty, 
we equip with a u;*-topology, induced by the seminorms T — > |Tr AT\ 
indexed by elements A e A. 

Lemma 3.2. Let 7„ G be a sequence of uniformly bounded density ma- 
trices, i.e. 

C := hmsup 1177111^1 = hmsup Tr SjnS < oo . 

n n 

Then one can extract a w* -convergent subsequence, 7^^. ^ 7, 7 G , and 
any w* -limit point 7 of the sequence {7n} satisfies ||7||-^i = Tr S'yS < C 

The convergence in Lemma 3.2 follows from standard application of the 
Banach-Alaouglu theorem. The positivity of the limit can be checked by 
testing with the projection operators Pf := |/)(/| G A for all f e H. 

Q.E.D. 



Define the following norms for operators 7^*^^ G C{H®^) 



7('^||„i,W :='IV 



Sx\ ■ ■ ■ Sxhl^ ^ ^Xk ■ ■ ■ S._ 



^Xl 



This norm is equivalent to the norm 



Tr|7(^) I + IVV.i v., . . . Va.,7^'=^V,., . . . V^^ 



and for density matrices we have the identity 

h^'^^ll^i.w =1V(/-A,J(7- A,J...(7- A,j7(') • 

Here we denote by W^'^*^) = rC{H®^) the set of operators with finite H}'^^'^ 
norm. Since higher derivatives on the same variable are not allowed, these 
norms are weaker than the operator analogue of the traditional higher order 
W^^^ Sobolev norms. 

The w*-topology on the space Ti}'^^^ is given by the seminorms indexed 
by the set 



A!^ ^ '■ — 'I'S'a;! • • • ^Xk ^ ^Xk ■ ■ ■ ^X\ 



K e 
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with norm 
and 



lU(fe) 



{Sxi ■ ■ ■ Sxf.) ^T{Sxi ■ ■ ■ ^xk) 
(fc) 



{•^^''\ II • lU(fe))* = (w^'^ II • ii7^i,{fc)^ 

analogously to the one variable case in Lemma 3.1. The analog of Lemma 
3.2 also holds. 



Define the set 

oo 

P:=|r = (7^,7(2),...) : jC^) en\H®'')] = Y[n\H^'') . 

k=l 

This set is equipped with the product topology r generated by the w*- 
topology on each li}{H®^). The convergence in D is characterized by the 
following property: r„ — > F in iV, r) as n — > 00, if for each k we have 
7^*^^ 7^'^) in w* sense in V}''^^''. 



Finally, we will consider time dependent density matrices. The set 
'(R+,(t)-Mt,W^'W) consists of functions 7W(t) : R+ ^ H^'C^) with 



supt(t) 



"i7^*'n*)llHi.W < 00. We define the set L°°(R+, = 
l^^ij^ (R+, {t)~^dt,'H^'^''^) and we equip it with the product of the w*- 
topologies on each factor. On each factor the w*-topo\ogy is given by semi- 
norms 

Jo 



T 



TV 



dt 



where A{t) G L\K+, (t)di,^(*^)), i.e. \\A{t)\\^ik){t)dt < 00. 



4 Structure of the proof 



In this section we give the main propositions and show how our theorem fol- 
lows from them. The proofs of the propositions are found in the subsequent 
sections. 



4.1 Cutoff of the Coulomb singularity. 



We now introduce an dependent cutoff for the Coulomb potential. De- 
compose V = Vi + V2 where 

Vi{x) = e{VNe-^x)Vix) 
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with a smooth cutoff function < 9 < 1 with 9 = 1 on outside of B{0, 2) and 
^ = inside .6(0, 1). Here B(x, r) denotes the ball of radius r centered at x. 
Notice the potentials, Vi and V2, depend on e and N which dependence is 
not explicitly labelled. 

We define the cutoff Hamiltonian 

1 ^ 1 

HN,e = - 2 XI + iV ^ ^^^^^ ~ 

i=\ e<j 
and the remaining part of the potential 

so that Hn = Hn,£ + W. 



The density matrices 7^^f ^ of the cutoff dynamics again satisfy the A^- 

(4.1) 



body Schrodinger hierarchy (2.5) if V is replaced with Vi, i.e. 
j^y{xi,...Xk;x^,...Xk) 

k rt 

^N,e,kithNo{xi, ...Xk] x[, . . . Xfc) + (-i) / ds UN,e,k{t " S 

1=1 -^0 

j dxk+i{yi{xi - Xk+l) - Vi{x'^ - Xfe+i)j 



N-k 



N 



£,(fe + l)/ / / X 

^^Ns [xi,--- Xk, Xk+i;xi, . . . Xfc, Xfe+ij 



where UN,s,k{t)l = ^•'^-fe and 



Wc now consider the evolution of the cutoff dynamics with smooth initial 
data. We introduce the notation 

N 



e=i 



We are ready to state the main propositions on the BBGKY hierarchy: 
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Proposition 4.1. [Apriori bound] Assume that the initial data of the N- 
body Schrodinger hierarchy (4-1) 

TiV.o := (tSJ^o' 7S,'i, . . . tS, 0' 0> • • • ) ^ C (4.2) 

forms a consistent sequence of normalized density matrices, i.e. = 
Trfc+i A; = 1, 2, . . . , A'' — 1 and Tr 7j^Q = 1. Assume that 

TrL^^^^^^<vl k = l,2,...N (4.3) 

for some constant i^o > independent of N. Let e > be any fixed constant 
and 

r?.,:=(7^J\7^!?,...) (4.4) 

be the solution to the N-body Schrodinger hierarchy (4-1) with initial data 
rjv,o oind Hamiltonian Hn,s- Then there exists a u > depending only on 
vq, and for any k fixed there is a constant N(e,k) such that r|^^ satisfies 
the a priori estimate 

||7^f^llHi.«<^' (4.5) 

for N > N{£,k). Therefore, Vj^^ forms a pre- compact sequence in L°^ (R.^,V) 
(with the w* topology) by Lemma 3.2. 

For ly > define the norm for T = (7^^) , 7^^) , . . . ) eV 

00 

k=l 

Proposition 4.2. [Convergence to the infinite hierarchy] Assume the initial 
data (4-2) satisfies the estimates (4.3). Lete > be any fixed constant. Then 
any vf limit point Too,t of the sequence \T%h\n=\,2,... C V satisfies the 
infinite BBGKY hierarchy (2.6) and for some v large enough ^oo,t\v is 
uniformly bounded for all t. 

At this stage we do not know that the hmit point Too,* = {l^oo,t^ 7oo!t! • • • ) 
forms a consistent family of density matrices. This will follow from the 
uniqueness of the solution to the nonlinear Schrodinger equation (2.10) and 
the uniqueness result below. 

Proposition 4.3. [Uniqueness of the infinite BBGKY hierarchy] Let T > 
be any fixed time. Then for any v > Q the infinite BBGKY hierarchy (2.6) 
has at most one solution in the set L°°([0,T], (P, || • ||;/)). 
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4.2 Approximation by smooth initial data 

For any k > let ipQ := e'^^ipQ. This smoothed initial data satisfies the 
following bound 

Uo-ro\\<CK\\M\H^ (4.6) 
using a simple estimate in Fourier space 

For S>0, let 

N 

-^Nfi{xi, ...xn):=YI ^o^i^j) 

denote the S/N regularized initial wave function. Let 7^^^^ be the projection 
matrix onto ^^q. The main reason to regularize the initial data is the 
following bound, whose proof is postponed. 

Proposition 4.4. There exists a constant C such that for any (5 > and 

k>0, 

= {^%,o,L'^%,o)<{c\\Mhy (4.7) 
whenever N is sufficiently large depending on S and k. 

The proof of this Proposition will be given in Sect 5. From now on we 
assume that the onc-particle wave function of the initial data is ^po G i?^(R^). 
We can now apply the previous setup of the e/s/N cutoff dynamics (Section 
4.1) to the regularized initial data. Let 

he the solutions to the Schrodingcr equation (2.1) with the original Hamilto- 

nian and with the regularized Hamiltonian, i^Af,e) respectively. Denote the 

5 (k) 6 £ (k) 

corresponding density matrices by 7;^ ^ and ■ From Propositions 4.1- 

4.4 and from the fact that F^^ (2-14) solves the infinite BBGKY hierarchy 
(2.6), we have the following corollary. Notice that the BBGKY hierarchy is 
independent of the cutoff e/\fN and the initial data for the BBGKY hier- 
archy is independent of the smoothing parameter 5/N . Both regulatizations 
disappear after we have taken the limit N ^ oo. 



Tr 



^k UN) 
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Corollary 4.5. Suppose that tpQ £ i7^(R^) and let tp^ be the solution to 
(2.10). There exists a constant v depending on \\iPo\\h^ such that for any k 
and 6 > 

hNf^\\nu>^)<'^'' (4.8) 
for sufficiently large N > N{S, e, k). Furthermore, 

N,t ^ i't 

as N ^ oo, in the w* -topology o/L°°(R_|_, (i)~^di, P). In other words, for 
each k 

S,e,{k) ^ Ak) 
'N,t ^ 'iPt 

in the weak* topology o/L°°(R+, {t)-^dt,V}{H®^)) as N ^ oo. 

We now control the deviation due to the cutoff and smoothing. The next 
Lemma shows that the A^-body wave functions for the 5/N regularized and 
the original initial data are close for all time. 

Lemma 4.6. Suppose the initial data tpQ G i7^(R^). Then we have 

limlimsupsupll^'Art - ^'ATtll = . (4.9) 

5-^0 N^oo t 



Proof. By unitarity, it is sufficient to consider t = 0; 



N 



N-j-1 



\\^N,o - *^,oll < E ll^oll'"' 11^0 - W II %^ < SWMrn 

using (4.6) and ||V'oll < IIV'oll = 1- Q.E.D. 



Finally, the following lemma controls the difference between the original 
dynamics and the £:/-\/]V-regularized dynamics. 

Lemma 4.7. For any e fixed, we have 

suplimsup II 5'^''^ - *Art|P < C{ipo)et , 

S N-*oo 

where C{iI)q) depends on the H'^-norm ofijjQ. 

The proof of this Lemma will be given in Section 8. Combining these 
Lemmas, we have 

lim lim limsupsup(t)~^||^'^^j — ^'Ar,t|| = . 

£^0 5^0 AT^oo t 
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Thus we have for any A; > 1 

lim lim limsupsup(t)~^||7f)-^/'^^ — 7l^l||i = . 

£-»0(5-*0 N-^oo t 

Since the trace norm is stronger than the w*-topology in C^{H^^'''>), The- 
orem 2.1 follows from Corollary 4.5 and Lemmas 4.6 and 4.7. Notice that 
the limit in Corollary 4.5 was controlled in the w* Sobolev norm. The con- 
trol in Lemmas 4.6 and 4.7 was in strong sense but without control on the 
derivatives. This is why the final result is only in the w*-topo\ogy of trace 
class operators. Q.E.D. 



5 Apriori bounds 

In this section we first prove Proposition 4.4 then Proposition 4.1. Recall 
that Proposition 4.4 guarantees that the smoothed initial data satisfies the 
apriori bound (4.3) for the initial condition in Proposition 4.1. 



Proof of Proposition 4-4- Let n = (ni, n2, . . . n^) be a sequence of positive 
integers, let k{n) = ni + . . . + ne and £{n) = £. We also define Lj := 
jf{I — ^xj), clearly L = Ylf=i ^j- "^^^ following identity is easily checked 

k N 

^' = E E E ^^■■■i^:- (5-1) 

£=1 n ; k{n)=k : i{n)=i ii,...i^=l : disjoint 

We have (recall k = S/N) 

(v'o^(^-A)>o") = 

< 
< 



jil+p^re-^P'\Mp)\'dp 
C^{l + n\K''-^)\\Ml2 . (5.2) 



18 DERIVATION OF NONLINEAR SCHRODINGER EQUATION . 



Therefore 

E {^%,o,Ll^LlK..Ll^^i^,) (5.3) 

il,...i£=l ; disjoint 

m ^ 



k ^ 





k 



k { I \i l = k (i.e. m = 712 = . . . = 1) 



A:! ((5-'=)iV-i if ^ < ifc 



Hence from this estimate, (5.1) and from the exponential bound on num- 
ber of possible n sequences with fc(n) = A; we obtain 



[^^,o,^'*Ar,o) < (1 + k\ {d-')N-^) (5.4) 

and (4.7) follows. This completes the proof of Proposition 4.4. Q.E.D. 



Proof of Proposition 4-1- We define 

1 ^ 

i=l 

where the constant /? > 2 will be chosen later. We also introduce Uij := 
N-^Vi{xi - Xj) (for i ^ j), Uii = and [/ := Ei<i<j<iv Uij. Thus 

HN,e + P^N = N{Q + U) . 

We also define 

N 

^ '•— ^ ^ QiiQi2 ■ ■ ■ Qik • 

U i*2v*fc=l • disjoint 

Note that Q,Q^^^ and U depend on N, but this fact is suppressed in the 
notation. The key Proposition is the following result to compare with 

trfe . 
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Proposition 5.1. There exist a constant (3 (depending only on fi) such that 
for any k and any N > N{s, k) we have the following operator inequalities 

C^Q'' <{Q + Uf < C^Q^ (5.6) 

with some positive constants Ci , C2 . 

First we show how this Proposition impUes Proposition 4.1. The assump- 
tions (4.3) and (5.6) imply that 

iov k = 1,2, . . . N ii N is large enough. Since any function of the Hamiltonian 
is conserved along the time evolution, we have 

TV 7^!r(^iV,. + P'^Nf < Cl{u^ + /?2'=)iV^ . 

By (5.6) and Q^^^ < we obtain 

^7^;fQ(')^C|(-o' + /3^')- 
Using the total symmetry of the TV-body density matrix and P > 1, 

and (4.5) follows. The pre-compactness of the sequence follows from 
the multivariable analog of Lemma 3.2 and a diagonal selection procedure. 

Q.E.D. 



Proof of Proposition 5.1: We proceed by a step two induction on k. The 
statement (5.6) for A; = is trivial. For k = 1 the bounds follow from 

\U\ < C(3-^Q (5.7) 

which is a consequence of Hardy's inequality in three dimensions (Uncer- 
tainty Principle Lemma in Vol.11. Section X.2. of [7]) 

^<-A., xeR'. (5.8) 

Assuming (5.6) is valid for some k, we show it for k + 2. Writing (Q + 
i7)'=+2 = (Q + U){Q + Uf{Q + U) we have 

(Q + U)Q^{Q + U)<{Q + Uf+'^ < CliQ + U)Q''{Q + U) . 

In order to compare UQ^U with Q^^'^ we will prove the following lemma: 
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Lemma 5.2. For < k < N/2, < a < 3 

UQ'^U < (^CP-^ + C{k) [/32feiV-'= + N^-a/2^-2k-2+a ^ ^-ij ^ gfc+2 _ 

(5.9) 

with constants depending on a and we can set C(0) = 0. 

Remark: A similar estimate holds for potentials much more singular than 
the Coulomb potential. With essentially the same proof, one can establish 
it for V{x) ^ \x\~^~^^ for any 77 > 0. 

From the Schwarz inequality A*B + B*A < A*A + B*B, with A* = 
VlUQ''/^, and B = -^Qi'/^+\ we have 

- (UQ''Q + QQ^U^ < 2UQ^U + ]^QQ^Q . 

Together with Lemma 5.2 with some 2 < a < 3 and (3 large enough, the 
lower bound in (5.6) follows from 

(Q + U)Q^{Q + U)> ^QQ'^Q - UQ^U > ^(5*^+^ 
asN>N{e,k). So Ci = 1/2. 

For the upper bound, we again let 2 < a < 3. Similar arguments lead to 
(Q + C/)Q*^(Q + U)< 2Q^+^ + 2UQ^U < CQ^+^ 
for N large enough, which completes the proof of Lemma 5.1 Q.E.D. 

Proof of Lemma 5.2. From Schwarz inequality, we obtain 

N 

For k = 0, (5.9) follows from 

Pijf < CN-^Qj < Cp-'^N-^QiQj , (5.11) 
which is a consequence of Hardy's inequality (5.8). 
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For k > 1 and for each fixed i,j we write Q = Qi + Qj + Qfj. By the 
weighted Minkowski inequality, we have for some constant C 



Q''<{Ck)\Ql + Q';) + 2Ql-.. 



(5.12) 



Let Di := \/—lN ^^'^Vxi- We have from Leibniz rule and Schwarz in- 
equality that there is a constant C{k) such that 



\Uij\ i^f \Ui^\ < CDf \Uij\'Df + C{k) L>f— {Dr\Uij\) 



m=l 



To estimate the derivative of Uij we notice that 



k—m 



2m+2-a 



(5.13) 



(5.14) 



for a < 2m + 2. We will choose 2 < a < 3. The singularity of |x — y| " is 
controlled by the following lemma whose proof is postponed: 

Lemma 5.3. Let W e ^^(R^) be any nonnegative potential, then 

W{x -y)< C\\W\\li{I - A^)(I - Ay) (5.15) 

on L2(r3 X R3). 



Combining this lemma with (5.14) we obtain (m > 1) 



(5.16) 



We have < C(A:)[L>f + (/J^Af-i)^]. Hence from (5.11), (5.13), and 
(5.16) we obtain 

\Uij\ \Uij\ < c(A;)iv-2(^(/32iv-i)'=+ivi-«/2£-2fc-2+«+iv-i)g,-(gf+g0 . 

Using g = Qj, Ei Qi < Q' and Q>I, 

N 



(5.17) 



22 DERIVATION OF NONLINEAR SCHRODINGER EQUATION . 



Finally, we use that Uij and Q^. commute and the estimate (5.11) to bound 



N 



N 



Yl Qfj ^^^Yl Cp-^N-'^QiQjQ'k < CI3-^Q^+^ . (5.18) 

i,j=l i,j=l 

Hence from (5.12), (5.17), (5.18) and Q > / the statement (5.9) follows. 
This finished the proof of Lemma 5.2 Q.E.D. 



Proof of Lemma 5.3. Let 



Hix) := 



dp 



(5.19) 



R3 a/1 

then H(x) ^ around the origin, smooth outside and it decays faster 

than any polynomial at infinity. Clearly (5.15) is equivalent to 

1 . 1 



V(/-A,)(I-A, 



--W{x-y) 



V(/-A,)(/-A,) 



< C||W||ii . (5.20) 



Let / G L2(R3 X R3), then 
1 



V(/-A,)(/-A,) 



Wix - y) 



^{I-A,)il-Ay) 



f 



(5.21) 



dxdx'dx"dydy'dy" f{x, y)H{x — x')H{y — y')W{x' — y') 
xH{x' -x")H{y' -y")f{x",y") 



dxdx' dx" dydy' dy" | / (x. 



\x — x \ \y — y \ 
\x' — x"\ \y' — y"| 



3/4 



H{x - x')H{y - y')W{x' - y')H{x' - x")H{y' - y") 

by symmetry of {x^y) and {x",y") and a Schwarz inequality. By the prop- 
erties of H 



|n|3/4 

hence x", y" integrals can be performed. We also define G{u) := |wp/^|i?(?x)| 
and notice that ||G||l2 < C. Finally, by Young's inequality: 

sup / G{x - x')\W{x' - y')\G{y' - y)dx'dy' < \\G\\12\\W\\li 

x,y J 

we obtain (5.20). Q.E.D. 
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6 Proof of the convergence 



In this section we prove Proposition 4.2. The uniform boundcdncss of 
||roo,t||i/ foUows from (4.5) and Lemma 3.2. To check that the hmit sat- 
isfies (2.6), it is sufficient test the equation against trace class operators 
O <E L^{H'^^) on a fixed time interval [0,T]. We note that satisfies the 
finite Schrodinger hierarchy with the cutoff potential (4.1). 

Since O is compact, we have 

Jim TVO(7-J)-7S)=0. 

This shows the convergence of the left side of (4.1) to that of (2.6). 

For the convergence of the first term on the right side of (4.1), we now 
show that UN,k,e{t) converges to Uk{t). By one-step Duhamel expansion we 
have 



j=i e<j 



(6.1) 



3=1 



From the triangle inequality, the unitarity of e*^* '^^^N,e ^nd \Vi \ < Ce ^N^^'^, 
we can bound the operator norm of the left side by 



(ne-^) 



3=1 



...Jo 



< 



e<3 
CPs 



Using (6.2) the convergence of 7^0, the fact that Tr7^''Q < Ttj'^'q = 1 and 
< 00, we thus have 



du 



(6.2) 



lim 



TtO 



(fc) 



. 



We now show the convergence of the second term on the right hand side 
of (4.1). The prefactor [N — k)/N can be replaced by 1 in the N ^ 00 limit. 
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since 



J Tr OUN,k,e{i-s) J dxk+i(Vi{xe-Xk+i)-Vi{x'^-Xk+i)^'y%\ 



(fc+i) 

s 



ds 



< ^IPIIIlFilloo sup IV 7^^+^) = 0(iV-V2) ^ 

se[o,t] 

for any fixed k and e using that ||Vi||oo = O^N^I"^). We then show that 



-Uk{t - s) 



J dxk+i (v{xi - Xk+i) - V{x'ii - Xfe+i)^ 7, 



(fe+i) 

CXD.S 



ds = 



for any i<k. Decompose this difference into three terms (/.) + (//.) + (HI-) 
where 



(6.3) 
ds , 



j dxk+i {Viixi - Xk+i) - Vi{x'^ - Xjk+i)^7^[^"'"^^ 
(77.) = j\oUk{t-s) J dxk+i(v{xe-Xk+i)-V{x'i-Xk+i)) 

ds , (6.4) 

(777.) = - [ TrOUk{t-s) (6.5) 
J dxk+i(v2ixi - Xk+i) - Vi{x\ - Xjk+i)j7^[^"'"^^ 
We estimate these three terms separately. 
The first term is bounded by 



-(its"'' -lit') 



ds 



|(7.)| < tsup OUj^^kS-s)-OUk{t-s) 

s<t 



X sup Tr 

s<t 



J dxk+i(Vi{xi - Xk+i) - Vi{x'^ - Xk+i)^-f%[ 



(fc+i) 

s 
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The norm \\OUN,k,£{t — s) — OUk{t — s)\\ is bounded by 



|CWjv,ifc,£(t-s)-OWjk(t-s)|| < ||C|| sup 

s<t 



~ eVN 

where we have Tiscd (6.2) in the last inequality. To bound the trace term, 
from the triangle inequality we reduce it to estimate 



Tr 



e,(k+l) 
N,s 



for s < t. Recall that we interpret the X|^_^_l integration as a partial trace 
over this variable. We can now use (3.2) and bound this term by 

Tr\Vi{xe-Xk+ih'£^^^\ 
where the trace is now over all + 1 variables. We define 

^ '.= Sx'^Sx2 ■ ■ ■ • 

From the inequality (5.8) we know that 

\V^{xt - xk+i)\^ < \Vixe - xfc+i)|2 < CSl^^ < C[S^^+y . (6.6) 
Hence for any self-adjoint 7 G we have 

Tr 



Vi{xi - Xk+ih = TI^J-iV^-i < CTr^7[5('=+i)]27 = CTr^5(^+i)725(^'+i) 

(6.7) 



using the cyclicity for the trace of the square root (9.1). Since the square 
root is monotonic in operator sense and 5'*^'^+^^ > /, we obtain 

CTry^5(^+V^S(^ < CTVy^ 5(fc+i)^[5(fe+i)]2^5(fe+i) = C\\-i\\^uk+i) , 

(6.8) 

where the last equality is the definition of the norm. Therefore we showed 
that 



Tr 



j dxk+iVr{xi - Xfc+i)7^J+'^| < C\\i 



N,s 11^1.(^+1) 



which is uniformly bounded by (4.5), hence the first term (I.) vanishes as 
N ^ 00. 
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The first half of the second term {II.) is written as 



where the trace refers to aU A: + 1 variables. In order to check that (//.) van- 
ishes in the limit ^ oo, it suffices to show that A{s) := [S^''+^^]-'^OUkit- 
s)V^,fc+i[5('=+i)]-i is in the space L°°([0, t], /C(''+^)), i.e. that it is compact 
and is uniformly bounded for all < s < t. Here V^^fe+i is the multiplication 
operator by V{x£ — x^+i). 

For uniform boundedness, we only have to show that Uk{t — s)V£^k+iS^^i 
is uniformly bounded, since O is bounded. Here Sk+i ■= Sxi^_^_i for brevity. 
Prom the unitarity of Uk{t — s) , it suffices to prove that 



This follows from the inequality (6.6). 

For the compactness, it is sufficient to show that S'j^^^OUk{t— s)V£^k+iS^^-^ 
is actually Hilbert-Schmidt for all s. Recall that O is a compact operator on 
H^'^, i.e., in the space of the first k variables, and the unitary map Uk{t — s) 
acts trivially on the (A;-|-l)-th variable. Hence we can define Og '■= OUk{t—s), 
which is a Hilbert-Schmidt operator on H^^. To emphasize that OUk{t — s) 
is viewed as an operator on we will write Og Ik+i- The identity 

Ik+i acts on the (k + l)-th variable. 

We introduce X := {xi, . . . Xk) and Y = (yi, . . . Uk), and we compute the 
Hilbert-Schmidt norm of S^}^^OUk{t — s)Vi^k+iS^}^i- 




< C . 



(6.9) 



Tr{Os ® Ik+i)S^l^Vi,k+iS^l^Vi 

= Tr(Ps ® Ik+i)Vi,k+iS^l^Vi,k+iSkli{Ps ® Ik+i) 




1 



{Vk+i - Xk+i)Os{Y; X)dXdYdxk+i^yk+i 



X 



/-A 




2 



1 



{Vk+i-Xk+i) dXdFdxfc+idyfe+i 



X 



/-A 




(6.10) 



L. ERDOS, H.-T. YAU 



27 



since Os is Hilbert-Schmidt. We used a Schwarz inequality in the third Hne 
and the fact that 

where (/— A)^^(.t— y) denotes the translation invariant kernel of the operator 
(/ — A)^^. This completes the estimate of the second term (//.)• 



dz 



ap 

n + » 



|2^2 



< oo 



Finally the third term is estimated as 

(///.) < 2t\\0\\ sup TV / dxk+iV2ixi - Xk+ih 

s<t J 

From (3.2) and Lemma 5.3 we can bound the trace term by 



£,(fe+l) 
N,s 



Tr 



where we introduced 7 := 'J^^g'^^^ for brevity and we computed ||V"2^||li < 
CeiV-1/2. Clearly 



T^\f^SfSl^ < IV^7[5(^+i)]27 < Chilli. 



(fc+i) 



as in (6.7)-(6.8), which is uniformly bounded in N by (4.5). Therefore {III-) 
goes to zero as N ^ 00, and we have proved Proposition 4.2. Q.E.D. . 



7 Uniqueness of the infinite hierarchy 

Proof of Lemma 4-3- Since the equation is linear, it is sufficient to show 
that the solution F^ up to time T is identically zero if the initial condition 
is Fq = and if for some K the apriori bound ||Fi||i, < K holds for t < T. 
Notice that it is sufficient to show that uniqueness holds for a short time, 
t < T{u), then one can extend it up to time T since the apriori bound holds 
uniformly for t <T. 

We need the following lemma whose proof is postponed. 
Lemma 7.1. For any 7 G 

k>l ande<k 



Trx 



J dyy(x,-2/)7(X,y;X',y)5('=)| <C||7||wi.(.+i) • 



where X = {xi, . . . x^) for brevity and recall S^'^^ = S'xi«S'^2 ■ ■ ■ 
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Armed with this estimate, we obtain from (2.6) and 'y^^'' = that 

hfLu.) = TV I^W^f 5(^^)1 (7.1) 

< 2^ Td^Tr 5«TV,,^Jy(a;,-Xfe+i)7f+^)]5('^) 
i=i 

,(fc+i) ■ 

Jo 

Let A{k,t) := Ht*'^^ ||-^i,(fc) • After iteration, we have 

ft fSl fSn-l 

A{k,t) <C"'k{k + l)...{k + n-l) dsi ds2 • • • / ds„ A(n + A;, s„) . 

Jo Jo Jo 

(7.2) 

Since Tt is bounded in (P, || • ||i,), there exists some K such that 

A{n, s) < Kv"" 

for all < s < T. After using this estimate in (7.2), the multiple time 
integration is bounded by f^/nl. This n! compensates the combinatorial 
prefactor k{k + 1) . . . (fe + n — 1) and the result is 

A{k,t) < if(2z/)"+'=(C0" < K{2iy)''{Cut)'' 

for any n. If t < T{u) := (Cz^)~^, then letting n — oo we obtain that 
A{k,t)=0. Q.E.D. 



Proof of Lemma 7. 1 . Without loss of generality we can assume that 1 = 1 
and we set x = xi, i.e. X = {x,X2, ■ ■ ■ x^). We define S := 8x28x3 ■ ■ ■ Sx^ 
and let 7 := SjS. Let V be the multiplication operator by the function 
{x,y) ^ V{x - y) on L2(Rf x R^) and let = {Px,i,Px,2,Px.^) be the 



components of momentum operator in the x variable, i.e., P^j := 
is similarly defined. 



-iVi. P, 



We denote R := Try{Vj) the partial trace of Vj, then we have to esti- 
mate 



Tlx 



S« J dy y(xi-y)7(X,y;X',y)S« 



= Trx 



Sx-RSx 



Using the identity = I+Yl^j=i ^xj Lemma 9.2 from the Appendix, 
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we estimate 



< CTrx 



SxRS^RSx 



SxR^Sx 



1/2 
1/2 



(7.3) 



SxRPxjRSx 



1/2 



CTi:x\RSx\ + CY,T^x 



Px,j RSx 



Hardy's inequality (5.8) implies that 

V"^ <Sl, < ^2 



(7.4) 



In the first term of (7.3) we use (3.2), (7.4) and the fact that Sx and Sy 
commute to obtain 



< Tr 



Ttl^SxW^S, 

X 1/2 



1/2 



< Tr I^SxjSIjSx) 

< TV (^Sx^SlS^^Sx) = TV (^SxSyjSljSxSy) 

< Tll^SxSy'jSxSy'-fSxSyj = TllSxSyjSySxl 

= Il7llwi.(fc+i) 
where TV denotes trace in all variables X, y. 

For the second term in (7.3) we first notice the following identity 

Px,jV=-PyjV + VPx,j + VPy,J 

Using TV|A + B\< Tr\A\ + TV|B|, we obtain for each j = 1, 2, 3 



TVx 



TtyPxjVjS, 



< Trx 



(7.5) 



+TVx 



TryVPxjjSx 



+ TVx 



TryVPy^jjSx 



< TV 



yjSxPyj 



+ TV 



VPxjlSx 



+ Tr 



VPyjlSx 



Notice that in the first term we used the cyclicity of the partial trace (9.4) 
from the Appendix before applying (3.2). 
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The estimate of the three terms in (7.5) are similar by using the inequal- 
ities (7.4) appropriately. The estimate of the first term is 



Tr 



1/2 



1/2 



(7.6) 



where we used the cyclicity (9.1) several times. 
The second term in (7.5) is estimated as 



Tr 



1/2 



= Tr [SxjPo:,jV^Pa,,jjS^^ 

< Tr (^SxjSlS^y^S^y^' = Tr (^S^Sy^Sl^S^Syf^^ 

< Tr(^S,Sy^SlS'y^S,Syy^' 



= llTllwiXfc+i) • 

The estimate of the third term in (7.5) is identical just the second inequality 
(7.4) is used. This completes the proof of Lemma 7.1. Q.E.D. 



8 Removing the Coulomb cutoffs 

Proof of Lemma 4-7: We fix N,5,e and for simplicity let := ^f^^, '■= 
'^^j^^. Clearly = $o- We compute 

dt\\^t-%f = 2Re {^t-^t,wyt) 
< 2||*t - ^'tll ||Ty*t|| . 
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We need to estimate 

\\W^tf (8.1) 

= ~ J \^2{xi-X2)\'^\^{xi,...XN)\'^dXi...dXN 

1 N{N -1){N -2) 
^iV2 2 

X J \V2ix1- X2)\\V2ix2- X3)\\^{xi,...XN)\'^dXi...dXN 

1 N{N-1) {N-2){N -3) 
^iV2 2 2 

X J \V2{xi- X3)\\V2{X2- X4)\\^ixi,...XN)fdXi...dXN 
< J \V2{xi-X2)\'^\^ixi,...XN)fdxi...dXN 

^'nJ\V2{x2 - 0:3)1 (|V^i*(xi, . . . Xjv)P+ l^ixi, . . . XN)\'^)dxi ...dXN 
J \V2ix1 - X3)\\V2{X2 - X4)\\^{XI, . . . XN)fdxi ...dXN , 

where we used that |V2(a;i — X2)\ < \xi — X2\~^ < C(/ — A^-^) for any X2 
fixed in the second term. 



It follows from Lemma 5.3 that 

X{\x -y\^< C{i,)X^-{I - A,){I - Ay) (8.2) 



\x - y[ 

as operators on L^(R'^ x R^) for A < 1, < k < 3. We use this estimate 
with K = 1,2, and A = ej^fN to continue the estimate as 



\W^t\\' < C-=l ||Va.iV^2^f + 11*1 



(4) 



for large enough A^. Using (4.8) this completes the proof of Lemma 4.7. 

Q.E.D. 
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9 Appendix on partial traces and density matrices 



In this Appendix we collect a few facts about trace class operators. We start 
with the analogue of the monotone and dominated convergence theorems. 

Lemma 9.1. Let < An < A be self-adjoint operators defined on a common 
separable Hilbert space. 

(i) If the sequence is monotone, An < An+i, An ^ A strongly, A is 
compact and sup„ Tr An < oo, then Tr A = sup„ Tr An, in particular A is 
trace class. 



(a) If An — ^ strongly and Tr A < oo, then Tr An — ^ 0. 



Proof: Let {fj} be an eigenbasis for A. For part (i), we have for any M that 

M M 

yZ{fi,Afi) = lim yZifuAnfi) < sup TV An . 

Taking the hmit M — oo we see that Tr ^ < sup„Tr An- The other 
direction is trivial. 

For part (ii), we notice that An is trace class and 

oo Af oo 

Tr An = Y,{fi, An fi)<Y,{fi, An fi)+ Yl ifuAfi). 
i=l i=l i=N+l 

The second term is smaller than any given ry > if A'^ is big enough since 
Tr ^ < oo, while the first term goes to zero for any fixed iV as n — oo. 

Q.E.D. 



Next we show the following properties of the trace of the square root. 

Lemma 9.2. Let A, B be positive self-adjoint operators on a common Hilbert 
space. Then 

Tr VAB'^A = Tr VBA^B (9.1) 

and 

Tr^/A + B <2(Tr^/A-\- Tr^/B) . (9.2) 
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Proof. We can assume that at least one side of (9.1), say TrV AB'^A is finite. 
Then X = AB^A > is trace class and so is Y := BA^B by cyclicity of the 
trace. By functional calculus 

/•oo 

VX = (const.) / t-^/^Xe-^^dt . 
Jo 

Since X is bounded, the operator Xe^^^ can be expanded into convergent 
power series for any t > 0. Using the cyclicity of the trace 

TV X** = Tr (AB^A)"" = TV (BA^B)"" = TV 

for any n, therefore 

TV Xe"*^ = TV Fe"*^ 

for any t > and (9.1) follows after integration. We used Lemma 9.1 (i) to 
interchange the trace and the improper dt integration. 



For the proof of (9.2) we let \i(A) > X2{A) > . . . > be the eigenvalues 
of A, similar notation is used for B and A + B. Fan's inequality [8] asserts 
that 

\n+m+M + B)< \n+l(A) + \m+l(B) 

for any m, n > 0. Therefore 

\2n+l(A + B) < \n+l(A) + Xn+l(B) 

X2n{A + B) < \2n-l{A + B)<Xr,(A)+Xn(B) 

and 

oo 

TVV^4+:b = J2Vh(A + B) 

k=l 

oo 

< 2j2VKiA) + Xn(B) 

n=l 

oo 

< 2J2{VK(A) + V\^) 

n=l 

= 2(TtVA + Ti Vb^ . 

Q.E.D. 



Now we give the definition of the partial trace. For more details, see e.g., 

[6]. 
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Definition 9.3. Let Hi,H2 be separable Hilbert spaces. Let A be a trace 
class operator on the tensor product H = Hi H2. Then there exists a 
unique operator B in C^{Hi) such that 

Tru, [BK] = Tru [A{K ® I)] (9.3) 

for every K G /C(ifi). This operator B is called the partial trace of A with 
respect to H2 and is denoted by TrH2-^ or just by Tr2 A. 

The existence of the partial trace follows from duality: the linear func- 
tional K TrHA{K (g) I) defines a bounded linear map on IC{Hi) and 
[)C{Hi)]* =C\Hi). 

Proposition 9.4. The partial trace satisfies the following relations 

Trn, I Trn, {I (S) A)B\ = Trn, \ TrH,B{I A) \ (9.4) 

and 

TrH,\TrH2A\ < TrH\A\ . (9.5) 

Proof. For the proof of (9.4) we use the variational principle (2.12) 
Tr 1 1 Tr2 (/ A) B I = sup Tr 1 ilT (Trs (/ A) 

= supTri,2 {K ®I){I® A)B 

= sup Tri,2 {K^I)B{I^ A) 

= supTri K(Tr2B(7(g)^)) 

= TVi|Tr2B(7®^)| (9.6) 

where the supremum is over all K G /C(i7i) with \\K\\ = 1. We used the 
cyclicity of Tri^2 and that / (g) A commutes with K ® I. 



For the proof of (9.5) we first observe that the variational principle (2.12) 
extends to bounded operators as follows: 



\A\\ 



sup 

LeC{H) : \\L\\ 



Ti AL 



(9.7) 



for any A G C^{H). The proof follows from /C(i7) C J0.{H) on one hand, and 
from \TrAL\ < ||A||i||L|| on the other hand, using (3.1). 



Therefore 



Tri |Tr2 A| = sup Tri L [Trg A] 

LeC{Hi) : ||L||=1 

= sup Tri,2 {L (g) I)A 

LgC(Hi) : \\L\\=1 

< TVi,2|A|. 
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The estimate follows again from the variational principle (9.7) since L®/ G 
C{Hi X H2) and \\L I\\ = \\L\\. Q.E.D. 



Finally, we prove that the space A is indeed the dual space of the Sobolev 
space H^: 

Proof of Lemma 3.1. (i) C A* . Let T G and we write any element 
AeAssA = SKS with some S e IC. Then 



TrTA 



Tr TSKS 



< \\K\\ Tr \STS\ = \\SKS\\a ||T||„i 



which shows that : ^ ^ Tr TA is a continuous linear functional on A 
indexed byT eH^ and |Kt|U* < \\T\\n^. 

(ii) .4* C H^. Let / G .4* be a continuous linear functional on A. 

Consider ip,(l) €^ i7^(R^) C H, i.e., Sip,S4> G H. Since is compact, 

\S(f)){Silj\ = S\(p){ilj\S is contained in A. Consider the sesquilinear functional 

J'f■.{1l;,c|>)^f[s\<^>){^P\S 

from X — > C. This sesquilinear map is continuous since 



/ s\4>){i^\s 



<||/IU* s\ci>){^p\s 



A* 



\4>){i^\ < 11/11 



\A* 



moreover it clearly extends from x to H x H. 

Hence there exists a unique bounded operator B = Bf such that 

{',p,Bcl>)=f[s\4>){'<P\s' 

for every (p,ilj & H, and ||i?|| < ||/||^* (uniqueness follows from the density 
of C H). 

By polar decomposition we write B = U\B\ with some partial isometry 
U. Let {ipi} be a finite orthonormal set, we can write 

N 



i 



i=l 



(9. 
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Here we used that U ipi is also orthonormal, and for any h E H 

i i 

i 

< \\hf. 

Hence taking the supremum for all orthonormal sets in (9.8), we see that 

i^l^l<ll/IU*- 

Finally, we define Tf := S~^BfS-'^ for any f eA*. Then clearly 

\\Tf\\n^=T^\Bf\<\\f\\A* 

and for any A e A 

Tr (TfA) = f{A) 

since this is valid for all A = 5"^"^-^ \ (/)i){ilJi\S finite range operators and 
these are dense in (A, \\ ■ m). Hence for any / G ^* we found a representative 
Tf G with smaller or equal norm. Q.E.D. 
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